Quantum Resonances of Kicked Rotor and SU (q) group 
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The quantum kicked rotor (QKR) map is embedded into a continuous unitary transformation 
generated by a time-independent quasi-Hamiltonian. In some vicinity of a quantum resonance 
of order q, we relate the problem to the regular motion along a circle in a (q 2 — l)-component 
inhomogeneous "magnetic" field of a quantum particle with q intrinsic degrees of freedom described 
by the SU(q) group. This motion is in parallel with the classical phase oscillations near a non-linear 
resonance. 
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For years, Chirikov's standard quantum map [[lj, i.e. a 
planar quantum rotor driven by very short periodic kicks, 
served as a cornerstone model for many investigations of 
onset and manifestations of dynamical chaos in quantum 
systems. Indeed, the standard map provides a local de- 
scription for a large class of dynamical systems and it 
is of quite general interest to understand its properties. 
The unitary Floquet transformation U which evolves the 
wave function ip(9) over each kick period is given by: 

U = U r -U k = exp (-\Trri 2 ) ■ exp (-ikcos6) (1) 

and consists of the operator Uk of a kick with strength 
k and a consecutive free rotation U r during the time T. 
Here rh — —id/dd and we put h = 1. 

In the classical limit the motion depends on the sin- 
gle parameter K = Tk. When this parameter exceeds 
unity the classical dynamics yields unrestricted diffusive 
growth of the mean kinetic energy. The quantum dy- 
namics is appreciably richer ||. In particular the phe- 
nomenon of dynamical localization has been discovered 
jlj , namely the quantum motion mimics the classical dif- 
fusive behaviour only up to a time tf> ~ k 2 after which 
it enters a stationary oscillatory regime. 

However, at present, we are still very far from a satis- 
factory understanding of this simple model of quantum 
chaos. The main obstacle for the analytical and numer- 
ical investigations are the so called quantum resonances 
which take place at a fixed value of k for the everywhere 
dense set of the rational values <; = T/Att = p/q, the inte- 
ger numbers p and q being mutually prime H,E| . In such a 
regime, the mean kinetic energy increases quadratically 
for asymptotic large times. Quantum resonances may 
strongly influence the motion in some vicinities around 
them. In other words they have certain widths. It is in 
this feature the main difference between the QKR and 
disordered quantum systems and the main reason for the 
failure of analytical tools in particular those which use 
supersymmetric technique ^J,^|. Therefore any attempt 
to improve our understanding of quantum chaotic motion 



must confront with the so far unsolved core problem of 
the motion in vicinities of resonances. 

In this paper we develop a general approach to the 
problem of the motion in a vicinity of a quantum reso- 
nance with arbitrary order q. We show that the most 
important role is played by the resonances with small 
q: for finite regions around them the quantum motion 
is explicitly shown to be regular and dominates the mo- 
tion for all ? values inside these regions. More precisely, 
the motion is described, for large though finite times, by 
a conservative Hamiltonian with one rotational degree 
of freedom and with a discrete spectrum. On the con- 
trary, when the resonance order is large enough the res- 
onant quadratic growth appears only in the remote time 
asymptotic and during a long time the motion reveals 
universal features characteristic of the localized quantum 
chaos. We demonstrate on a specific example that our 
approach allows to predict the quantum evolution with 
great accuracy for long time (see Fig.l) even though the 
corresponding classical motion is chaotic and exponen- 
tially unstable. 

Basically, we use the same idea of embedding the map 
(0) into a continuous unitary transformation which has 
been used before || in the study of the classical counter- 
part of (|l|). The Floquet operator is represented in the 
form U = exp(— iH), so that the motion appears as a 
conservative evolution in a continuous time t of a system 
with one rotational degree of freedom described by the ef- 
fective time- independent quasi-Hamiltonian Tl. The lat- 
ter cannot generally be found in a closed form. Rather, 
it is formally expressed as an infinite sum of successive 
commutators. However, when the condition of a quan- 
tum resonance is fulfilled it simplifies enormously (^] and 
reduces to a pure phase (or local gauge) transformation 
from the unitary unimodular group SU(q). Moreover, in 
some vicinities of quantum resonances, the sum of com- 
mutators can be truncated thus allowing to obtain new 
insight on the nature of the motion. 

Before turning to the consideration of the general 
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case, we illustrate our approach on the simplest and 
strongest resonances q — 1,2,4 (see also ||). For the 
main resonance q = 1 the rotation operator is equiva- 
lent to identity. The Floquet transformation ([!]) is then 
simply e~ m ( e \ where v(9) = kcos9 is the interaction 
potential. After t repeatitions of this transformation, 
an initially isotropic wave function gets ~ kt harmon- 
ics. This number is naturally measured by the operator 
rh(t) = e mt rhe~ lvt = m — v'(9)t. This yields the increase 
law Ek(t) = (Am 2 (i))/2 = rt 2 with the resonant growth 
rate r(k) = (v' 2 ) = fc 2 /4 for the kinetic energy. The 
symbol (...) stands here for the average over the angle 9. 

For the other two resonances q = 2, 4 the rotation op- 
erator U r has only two eigenvalues (1,-1) and (1, — i) 
respectively. The corresponding eigenfunctions tp^(0) 
are periodic and contain only even (odd) harmonics for 
the first (second) eigenvalue in each case. By writing the 
wave function in the form of a spinor ^f(0) 



E k (t)=rt 2 + X (t) 



(5) 



(transposed) ^ = ^ 



(2) 



the rotation is described by a matrix operator U r — 
e -ia r (a 3 -l) with ar = n / 2 for q = 2 and a r = -tt/4 
if q = 4. The kick matrix has in both cases the same 
form e~ mai . Here Uj are 2x2 Pauli matrices. Simple 
manipulations cast, up to an irrelevant constant phase, 
the Floquet operator, that acts on the spinor into 
the form of a transformation from the SU(2) group 

£/■(«») = exp (-iw n • er) = exp f-iH^ (6>)) , (3) 

where for q = 2, w — 7r/2, and n = (0, sin?;, cos v) 
while for q = 4, w = arccos(cosv/\/2) and n = (1 + 
sin 2 v ) _1//2 (sin?;, — sini> , cost* ). The hermitian operator 
in the exponent is equivalent to the Hamiltonian of the 
spin 1/2 in a magnetic field oriented along the vector n. 

The unified representation ^ allows us to calculate 
in closed form the evolution of the angular momentum 
operator M = diag(rh, m) 



AM(t) = e 



iw n (T ti\/\r„ — iw n (7 t 



Jo 

. sin 2wt , 1 — cos 2wt . 
ujinH — n H n x n | ■ cr. 



(4) 



The 3 matrices S a , a = 1,2,3 are generators of the ad- 
joint representation of the SU(2) group while the vectors 
n, n' and n' x n define an orthogonal basis in this space. 
At a given angle, the evolution (||) consists of two dif- 
ferent contributions. The linearly growing part appears 
because of inhomogeneity of the magnetic field while the 
periodic part is due to the spin rotation in this field. The 
corresponding kinetic energy evolution reads 



where 



r(T,k) = l((w') 2 ), X (q;k;t) = h n' 2 sin 2 wt). (6) 



When q = 2, the quantity w = tt/2 is independent 
of the angle 9 so that the quadratically growing term 
disappears and the function x(t) is purely periodic 



E k (t) = ±fc 2 sin 2 (±7rf) 



(7) 



Only the spin rotation contributes and the energy jumps 
between values and k 2 /4 when the time t runs over 
integers. On the other hand, for q = 4 both contributions 
exist. The function %(i) fluctuates with time approaching 
as const j\ft the finite positive value 



k 2 r 

— / d8 sin 2 8 [1 + sin 2 w]' 



(8) 



Let us now detune slightly from the exact resonance, 
T = 7 1 ( res ) + k. The Floquet operator looks then as 
U M («) = exp (-|kM 2 ) V p re q s) = exp (-{H{n)). Repre- 
senting the quasi-Hamiltonian as 



H = nii { - res \9) + K 2 Q(K) 1 
we come to the condition 



exp(-|KM 2 ) = 

T*exp\-iK [ dre^ iTi< '" )r g( K )e^ < '" )r 



(9) 



(10) 



The symbol T* indicates the antichronological ordering. 
One can formally solve this equation by expanding the 
operator Q(k) over the small detuning k. The effective 
quasi-Hamiltonian resulting from such a procedure ap- 
pears in the form of the series 

H = Kh& es \6) + X - {J, F 1 (9)} + + X - JF 2 (9)J + ... (11) 

which is, in particular, an expansion in powers of the an- 
gular momentum J = kM in which the zero-order term 
generates the resonant Floquet transformation. All op- 
erators in (|TT| ) are one or two dimensional matrices de- 
pending on the order q = 1 or 2, 4. Being a factor in 
front of the angle derivative, the detuning k plays in jll| ) 
the role of the dimensionless Planck's constant. Keep- 
ing only the terms written in eq. (^) explicitly, the 
quasi-Hamiltonian 7i, when q — 2,4, is formally equiva- 
lent to the Hamiltonian of a "particle" with the spin 1/2 
which moves along a circle in an inhomogeneous mag- 
netic field. The term linear in the angular momentum J 
mimics a sort of spin-orbital interaction. Due to the pe- 
riodic boundary condition, the eigenvalue spectrum {e} 
of H is discrete so that the detuning restricts the growth 
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of the energy to a finite maximal value. This justifies the 
expansion over the angular momentum. 

Comparing the results predicted by the quasi- 
Hamiltonian with exact numerical simulations of the 
original quantum map (Q), we found that already the ap- 
proximation (|ll]) describes quite satisfactorily not only 
the cut-off of the initial resonant growth (when it exists 
in the resonance regime) and the mean value of the ki- 
netic energy after the growth has been saturated but also 
delicate features of very irregular quantum fluctuations 
in the plateau region. The term quadratic in the angular 
momentum J turns out to be of the principal importance. 

As a typical example we show in Fig.l the results for 
the resonance q = 4. The points are numerical simu- 
lations with the exact Floquet operator (|l|) (only each 
4th point is kept in the main part and 500th in the in- 
set) while the solid line corresponds to the evolution de- 
scribed by the time-independent quasi-Hamiltonian (11). 
The condition under which the influence of the lowest 
omitted correction is weak gives an estimate of the width 
Ak of a resonance [Q. This leads to Ak ~ 1/k for the 
resonances q = 1,2. The resonance q = 4 is much nar- 
rower and Ak tx 1/k 2 . Both estimates are in agreement 
with the numerical data. Outside these intervals the ex- 
pansion becomes transparently divergent. In Fig. 2 the 
transient region of the motion near the resonance q = 4 
as it comes from the exact numerical simulations is pre- 
sented. Two quite different regimes are clearly seen. The 
crossover from the regular (lower plateau which is well 
described in terms of the quasi-Hamiltonian) to chaotic 
(upper line) regimes takes place in a rather narrow in- 
terval of the detuning. To explore the regularity domain 
and the adjacent region, we have fitted the height E p i of 
the plateau in Fig. 2 as a function of the detuning k, (see 
Fig. 3). The plateau height scales as k _1 in the regularity 
domain (k < 10~ 4 ), while in the "quantum chaos" region 
(k > 10 -3 ), the plateau height is scattered around the ex- 
pected localization value oc I 2 . Note that the resonance 
essentially suppresses the diffusion in the intermediate 
region as well. Higher corrections to eq. ( |i"T| ) may be 
relevant in this domain. 

Rigorously speaking, an expansion of the kind (|TT| ) 
cannot, in spite of the satisfactory agreement with the 
exact numerical solution, converge. Indeed, an infinite 
number of quantum resonances of large order hit the do- 
main Ak of influence of a strong resonance with a small 
order. In such a resonance of high order q, the Floquet 
operator is a SU(q) transformation 



U^ s) = cxp {-iw n • A) ; n 2 = 1 



(12) 



which is a plain extension of eq. (g). Here the matrices 
A a are the generators of the g-dimensional fundamental 
representation of the group, n • A = J2 a n a^a and n is 
a unit vector in the (q 2 — 1) -dimensional adjoint space. 
The transformation (O) depends on the angle 8 via the 



periodic functions w{&) and n(0) which in turn are ex- 
pressed in terms of the kick potential v{&). These func- 
tions satisfy a system of q 2 — 1 transcendental equations 
which, generally, cannot be solved analytically. Never- 
theless, some generic information can be obtained even 
without knowing the explicit solution. Indeed, the evolu- 
tion of the q x q angular momentum matrix is described, 
similarly to eq. (0) , by 



AM(t) = [M (0) i + M (1) (i) 



(13) 



The matrices M^ ) are expressed in terms of the q-1 zero 
modes of the matrix (n • A) which acts in the (q 2 — 1)- 
dimensional adjoint space, while the quasiperiodic (as 
long as the angle is kept fixed) matrix M"(i) is formed 
by the q(q — l)/2 pairs of mutually conjugate modes with 
finite eigenvalues ±£ a . These two matrices determine 
respectively the resonant growth rate 



r(p,?;*) = !£«»7«W] 



(0)l2\ 



and the asymptotic value 

<teb>0} 



Xoo = 2 £ >0> (C 2 In' • X {b) I 2 E Q Va \xi b) I 2 ) 



(14) 



(15) 



The quantities rj a in this formulae select the (2q-l)- 
dimensional active subspace of the adjoint space, which is 
reachable for the motion under the isotropic initial con- 
ditions. A small detuning from the considered resonance 
kills the unrestricted growth created by the resonance 
itself and gives rise to the q x q quasi-Hamiltonian ma- 
trix of the same structure as in ( |iT|) . The motion looks 
like that of a "particle" with q intrinsic degrees of free- 
doms described by the SU (q) group, along a circle in a 
(q 2 — l)-component "magnetic" field. 

Returning to the convergence problem, one sees that, 
independently of the number of the corrections taken into 
account, the quasi-Hamiltonian approach fails to repro- 
duce the unrestricted resonant growth in the resonant 
points within the domain of the influence of a strong 
resonance (i.e. with a small order) around which the 
expansion is performed. Indeed, the spectrum of the 
quasi-Hamiltonian is always discrete whereas the reso- 
nant growth implies continuous spectrum. 

However, as confirmed by numerical data, the growth 
rate r(p, q; k) is exponentially small when q noticeably 
exceeds the localization length ||. Therefore, the res- 
onant growth reveals itself only on a very remote time 
asymptotics. Qualitative arguments presented in Q con- 
nect this fact with the exponentially weak overlap of 
the neighboring localized parts of the globally delocal- 
ized quasienergy eigenfunctions. Exponential effects of 
such a kind, which are characteristic of the tunneling, 
are well known to be beyond the reach of perturbation 
expansions. That is why they cannot be described in 
the framework of the quasi-Hamiltonian method. The 
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latter reproduces well only those features of the motion, 
which are determined by the discrete component of the 
quasienergy spectrum, in particular, the function %(£) 
which attains its asymptotical value (15) much faster. 
Therefore inside the width of a strong resonance the time 
dependence of these functions is dictated for all weak ones 
by their strongest brother. 

On the other hand, if the order is very large, q » I, 
and the resonance lies in the region of typical irrationals 
being far from all strong resonances, already a very small 
detuning suffices for killing the quadratic growth with ex- 
ponentially small rate |L0| . At the same time, such a shift 
does not influence the function x(t) which reproduces on 
a large (though finite) time scale characteristic features 
of the "localized quantum chaos" Q . 
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FIG. 1. The energy E times k versus the number of kicks 
for q=4. Deviation of the theory from the exact map is indi- 
cated at the bottom. 




FIG. 2. The crossover region near the resonance q=A. 
Double log. scale is chosen. The resonant behavior on the 
initial stage is seen in both cases. 




FIG. 3. The plateau height E p i versus detuning k near the 
resonance q=4. The solid line corresponds to the theoretical 
relation E p i ex 1/k. 

The behavior becomes more complicated if a number 
of resonances with comparable and moderate orders are 
neighbouring and their domains overlap. The expan- 



sion near one of them forms a plateau which lasts un- 
til the quadratic growth in a next resonance of the same 
strength reveals itself so that the original expansion fails. 
However, the expansion near the new resonance cuts off 
this growth and forms a higher plateau until a next reso- 
nance comes to the action. Such a pattern of repeatedly 
reappearing regimes of resonant growth has been discov- 
ered in Q 

In conclusion, the concept of the time-independent 
quasi-Hamiltonian of a quantum map is proposed. The 
motion of the quantum kicked rotor in a quantum res- 
onance of order q is exactly described by a continuous 
transformation from the SU(q) group. The motion in 
some vicinity of this resonance is proved to be similar 
to that of a quantum particle with q intrinsic degrees 
of freedom along a circle in an inhomogeneous (q 2 — 1)- 
component "magnetic" field. The latter is in analogy 
with the classical phase oscillations near a non-linear res- 
onance. The resonances with sufficiently small q master 
the motion in finite domains near the resonance points 
though the widths of this domains rapidly diminish with 
the order. 
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